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Methods from soft-coUinear effective theory are used to perform the threshold resummation of 
Sudakov logarithms for the deep-inelastic structure function F2{x, Q^) in the endpoint region x ^ 1 
directly in momentum space. An explicit all-order formula is derived, which expresses the short- 
distance coefficient function C in the convolution _F2 = C 8 0g in terms of Wilson coefficients and 
anomalous dimensions defined in the effective theory. Contributions associated with the physical 
scales and Q^{1 — x) axe separated from non-perturbative hadronic physics in a transparent 
way. A crucial ingredient to the momentum-space resummation is the exact solution to the integro- 
differential evolution equation for the jet function, which is derived. The methods developed in this 
Letter can be applied to many other hard QCD processes. 
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I. INTRODUCTION II. FACTORIZATION FORMULA 



A generic problem in applications of perturbative QCD 
to collider physics or heavy-quark physics is to disen- 
tangle contributions associated with different momentum 
scales, and to resum large logarithms of ratios of such 
scales to all orders in the perturbative expansion. In 
processes containing hadronic jets, a scale hierarchy is 
created by the fact that the invariant mass of a coUi- 
mated jet is typically much smaller than the hard scale 
of the process (e.g., the center-of-mass energy). The in- 
tricate interplay of soft and coUinear emissions then leads 
to large Sudakov double logarithms. The resummation 
of these logarithms is conventionally performed in mo- 
ment space, and predictions for differential cross sections 
in momentum space are obtained by an inverse Mellin 
transformation. This procedure is cumbersome and of- 
ten leads to unphysical singularities, because the resum- 
mation formulae involve integrals over the Landau pole 
of the running coupling. These singularities are dealt 
with by means of ad hoc prescriptions, or by introducing 
artificial infrared cutoffs. 

In this Letter we develop an approach based on ef- 
fective field theory, which allows us to resum Sudakov 
logarithms for a large class of processes directly in mo- 
mentum space. The starting point is a factorization the- 
orem for the cross section, in which contributions from 
different momentum scales are separated in a transparent 
way. Evolution equations for the various components in 
the factorization formula are solved exactly in momen- 
tum space, in such a way that one never encounters inte- 
grals over the Landau pole. We illustrate the procedure 
with the example of deep-inelastic scattering. However, 
the same methods can be applied to many other hard 
QCD processes, such as Drell-Yan lepton-pair produc- 
tion, prompt photon production in hadron-hadron col- 
lisions, Higgs-boson production in gluon-gluon fusion, 
heavy-quark fragmentation, event shapes, and others. 
Technical details of our derivations are presented in [1]. 



We focus on the flavor non-singlet component of the 
structure function F2{x,Q^) in deep- inelastic scattering 
(DIS) of electrons off a nuclear target, e~ -I- N{p) 
e~ -\- X{P), denoting hy q = P — p the momentum of the 
virtual photon. We are interested in the region where the 
Bjorken scaling variable x = Q^/{2p ■ q) is near 1, so that 
there is a hierarchy of scales :s> Q'^{1 — x) 3> Aq^j-,. 
The intermediate scale Q^(l — x) « Af|- is set by the 
invariant mass Mx of the final-state jet. In this region 
the structure function can be written in the factorized 
form [2-4] (with /x/ the factorization scale and Cg the 
quark electric charge) 

Frix,Q')=J2el\CviQ^f,fr 

xQ' [ dCJ{Q\^^x),^if)rg%^,Iif)■ 

This formula is valid to all orders in perturbation theory 
and at leading power in (1 — a;) and t^QYil^'^'x- Here Cy 
is a hard matching coefficient, J is a jet function, and 
(\)^ is the non-singlet component of the quark distribu- 
tion function in the nucleon. As shown in [1], a simple 
derivation of the factorization formula can be given using 
the technology of soft-coUinear effective theory (SCET) 
[5] (see [6-8] for earher investigations). 

In SCET the hard function Cy is identified with the 
Wilson coefficient in the matching relation of the QCD 
vector current onto the unique leading-power current op- 
erator in the effective theory. To calculate the Wilson co- 
efficient one must compare perturbative expressions for 
the photon vertex function in the two theories. The cal- 
culation can be simplified by performing the matching 
on-shell, in which case all loop graphs in the effective the- 
ory are scaleless and vanish in dimensional rcgularization. 
The bare on-shell vertex function in QCD (called the on- 
shell quark form factor) has been studied at two-loop 



1 



order and beyond [9-11]. The form factor is infrared di- 
vergent and must be regularized. When the SCET graphs 
are subtracted from the QCD resuh, the infrared poles 
in 1/e get replaced by ultraviolet poles. To obtain the 
matching coefficient we introduce a renormalization fac- 
tor Zy, which absorbs these poles. At one-loop order 
this gives [6] 



Cv(Q^/i) = l + 



CpCts 



-X^ + 3L - 8 + — 



where L = ln((5^//L<^) and = otsip)- The two-loop 
expression for Cy can be found in [1]. The scale depen- 
dence of the Wilson coefficient is governed by the evolu- 
tion equation 



dCv{Q\^x) 
dlnjj, 



In — -|-7^(as) 



(1) 



where Fcusp is the universal cusp anomalous dimension 
of Wilson loops with light-like segments [12], which is 
associated with the appearance of Sudakov double loga- 
rithms. The quantity 7^ accounts for single-logarithmic 
evolution effects. The anomalous dimension can be ob- 
tained from the coefficient of the 1/e pole term in the 
renormalization factor Zy- Using the results of [11] it 
can be calculated at three-loop order [1]. 

The jet function J is defined in terms of the disconti- 
nuity of a vacuum correlator of two quark fields, made 
gauge invariant by the introduction of Wilson lines. It 
obeys the integro-differential evolution equation [13] 



dlnjjL 



ZTcuspK) ln^ + 27^(a,) 



2 r cusp ( 



Jo 



p2 _ pl2 



dp'- 



We encounter again the cusp anomalous dimension, and 
in addition a new function 7'', which has been calculated 
in [13] at two-loop order, and whose three-loop coefficient 
is determined in [1]. 



III. SOLUTIONS OF THE RENORMALIZATION 
GROUP EQUATIONS 

The exact solution to the evolution equation (1) is 

Cy((^^M) = cxp [2S{^lh,^J) - a^v{^h,^i)\ 

-i2 \ -o-r(Mh:l') 
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(2) 



where /x^i ~ Q is a hard matching scale, at which the 
value of the coefficient Cy is calculated using fixed-order 
perturbation theory. The Sudakov exponent S and the 
exponents a-y are given by 



Q( \ /■ J rcuspCa) f da' 

ar' 



^^reusp(a)^ (3) 



/3(a) ' 



and similarly for a^v, where f3{as) = dag /din 11 is 
the /3-function. The explicit perturbative expansions of 
these expressions valid at next-to-next-to-leading order 
(NNLO) in renormalization-group (RG) improved per- 
turbation theory are given in [1] . 

An important object in the derivation of the solution 
to the evolution equation for ,/ is the associated jet func- 
tion j, which has originally been defined in terms of an 
integral over the jet function followed by a certain re- 
placement rule [14]. More elegantly, it can be obtained 
by the Laplace transformation 

i(ln^,M)=^°°dp2e-^^^ J(p^M), 

where s = l/{e^'^Q'^). The inverse transformation is 

Using the evolution equation for the jet function we find 
that the associated jet function obeys 



dln/x 

- 2r 



,K) ln^+27''(a,) 



which is analogous to the evolution equation (1) for the 

hard function. Inserting the solution to this equation into 
the inverse transformation (4) we obtain 



J{P^,IJ) = exp [-4:S{iJ,i,ii) + 2a^j{ii^,i^)] 



(5) 



where 77 = 2ar{^J.i, fi), and 9^ denotes a derivative with 
respect to this quantity. The above form of the result is 
valid as long as > (i.e., fj, < jii). For negative ry the 
singularity at = must be regularized using a star 
distribution [1]. Relation (5) is one of the main results of 
this Letter. It relates J to the associated jet function j 
evaluated at a scale fn, where it can be computed using 
fixed-order perturbation theory. At one-loop order 



3{L,ii) = 1 + 



47r 



2^2 _ 3L + 7 . 



27r2 
3 



where in (5) the argument L is replaced by the deriva- 
tive operator dr^. The two- loop expression for j can be 
extracted from [13]. 
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IV. MOMENTUM-SPACE RESUMMATION 

We are now ready to write down a resummed expres- 
sion for the structure function valid to all 
orders in perturbation theory and at leading power in 
(1 — x) and Aq(jj-,/M^. When combining the results (2) 
and (5) the Sudakov exponents can be simplified. Intro- 
ducing the short-hand notation a^^ 
after a straightforward calculation 



a^j — a^v, we find 



F^%x, Q2) =^el \CviQ^ t^h)\' U{Q, /x/., f^^, /x/) 



V /X- /I 



7EV 



(6) 



where 



U{Q,fih,IJ'i,fJ'f) = exp[4S'(/i^,/Xi) - 2ajy{ij,h,iJ,i)] 

-)2 \ -2or(jU)i,(Ui) 



91 



exp [2a^4>{ni, lJ.f )] 



and as before rj = 2ar{f^i, fJ-j)- The remaining integral 
can be performed noting that, on general grounds, the 
behavior of the parton distribution function near the end- 
point can be parameterized as 



where > 0. This leads to the final expression 



J2„(^lx(l)'^^{x,iif) 



jq -Tq 
X 



In 



CviQ'^.lih)? U{Q, Hh, Hi, lif) 
Q\l~x) 



3^7Er,r(l + b(/X/)) 

r(l + 6(/x/)+r?) • 



(7) 



The exact all-order results (6) and (7) are independent 
of the scales ^ih and /Xj, at which the^ matching coefficient 
Cv and the associated jet function j are calculated. The 
answers simplify further if wc choose the "natural" val- 
ues Hh = Q and /Xj = — X (for fixed x). In practi- 
cal calculations the residual dependence on these scales 
introduced by the truncation of the pcrturbativc expan- 
sions of the various objects can be used as an estimator 
of yet unknown higher-order corrections. 

Above we have accomplished the resummation of 
threshold logarithms for F2 directly in momentum space. 
The resulting formulae are simpler than corresponding 
expressions in the literature (see e.g. [15]) in that they 
do not require a Mellin inversion and in that the de- 
pendence on X and Q is explicit. The right-hand sides 
of (6) and (7) can be evaluated at any desired order in 
resummed perturbation theory. Using currently avail- 
able results, it is possible to include terms at NNLO 



[1], which is equivalent to the so-called next-to-ncxt-to- 
next-to-leading double-logarithmic (N^LL) approxima- 
tion. The resummation is under perturbative control as 
long as (l — x) ^ Aqqp/Q^, since only then the interme- 
diate scale m ~ Qy/l — x is a short-distance scale. While 
the theoretical description thus breaks down very close 
to the endpoint, wc note that weighted integrals of the 
structure function over an interval xq < x < 1 can be 
calculated as long as Q\/l — xq is in the short-distance 
domain. 

It is instructive to compare our result (7) with the 

conventional approach to threshold resummation in DIS, 
which proceeds via moment space [2,3]. One defines 



F^M(^^)= tdxx''-^F^{x,Q' 
Jo 



) 



= C^(Q2,/x/)^e2<^-^+i(/x/), 

Q 

where the moments of (/)"^(^,/x) are defined in analogy 
with those of i^"^(a;, Q^)- For large values of N the inte- 
gral is dominated by the endpoint region {l — x) ~ 
The short-distance coefficient Cjv is decomposed as 

CN{Q'^,^J.f) = .9o((5^,M/) exp [Gat (Q^,^/)] , 

where the prefactor go collects all iV-independent terms, 
and the exponent is written in the form (see [15] for the 
most up-to-date discussion) 



G^f(g^M/) = / dz 
Jo 



1 ^N-l 1 



1 - Z 



(8) 



Aq{as{k)) + Bq{as{QVl^z)) 



The resummation for the momentum-space structure 
function ^2(2^, Q^) itself is obtained from that for the mo- 
ments F2.n{Q'^) by an inverse Mellin transformation. It 
is possible to show (see [1] for details) that the outcome of 
this procedure is equivalent, at any finite order in the per- 
turbative expansion, to the result (7) derived from effec- 



tive field theory, provided we identify Aq (a^ ) 
and 



- cusp 



K) 



^+12^ 



B,(a,) =7-^(a,)+V lnj(0,/x) 



12 3 



rcusp(^5) ; 



where V = d/dlnn^. It follows from this relation that 
the quantities Bg and 7^^ agree at first order in (as 
observed in [6]), but they differ starting from two-loop 
order. 

There are a few unpleasant features of the conventional 

approach which are worth pointing out. First, note that 
the integrals over the functions Aq and Bq in (8) run 
over the Landau pole of the running coupling as(/[x), 
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FIG. 1. Dependence of the resummed result for F2^{x,Q'^) 
on the hard (left) and intermediate (right) scales, at differ- 
ent orders in RG-improved perturbation theory: LO (dot- 
ted), NLO (dashed), and NNLO (soUd). We use Q = 30 GeV, 
X = 0.9, Hf =5 GeV, and &(/i/) = 4. 



introducing an infrared renormalon ambiguity of order 
■'^qcd/^'^x- such problem arises for the integrals (3) 
in our approach. The particular integral representation 
of the solution (8) results from the fact that in the con- 
ventional approach one solves a set of partial differen- 
tial equations derived by diagramatic methods instead of 
the RG evolution equations in SCET [2,3]. The use of 
RG methods for the resummation avoids the Landau-pole 
singularity in the exponent [16,17]. Secondly, when per- 
forming the inverse Mellin transform one needs to inte- 
grate the function Gat over along a contour parallel to 
the imaginary axis. This integration involves arbitrarily 
small physical scales [fc^J ~ Q'^/\N\, leading to a second 
encounter with the Landau pole. Different prescriptions 
to deal with this problem have been proposed in the lit- 
erature. In our approach integrals over the Landau pole 
never arise, because factorization and resummation are 
performed directly in momentum space. The singulari- 
ties appearing in the conventional approach are an arti- 
fact of the way the resummation of large logarithms is 
implemented and they cannot be used to assess whether 
a corresponding renormalon pole is present [17]. In our 
approach infrared renormalons are expected to affect the 
large-order perturbative behavior of the matching coeffi- 
cients C'v and j. The corresponding infrared ambiguities 
will be commensurate with power corrections from sub- 
leading operators in the effective theory. The evolution, 
on the other hand, is driven by anomalous dimensions, 
which are expected to be free of renormalons. 

Figure 1 shows the scale dependence of our result for 
F^^ obtained with Q = 30 GeV, x = 0.9, Hf = 5 GeV, 
and n/ = 5 light flavors. Varying the scales /dh and fii 
about their default values, we observe that the residual 
scale dependence is strongly reduced when going to suc- 
cessively higher orders in perturbation theory. In the 
literature the matching scales arc often held fixed, and 
the perturbative uncertainties can only be estimated by 



comparing results at different orders in the expansion. 



V. CONCLUSIONS 

Using methods from effective field theory we have in- 
troduced a new approach to the resummation of large Su- 
dakov logarithms in hard QCD processes. Factorization 
and resummation are performed directly in momentum 
space, such that the resulting formulae are free of un- 
physical infrared sensitivities and provide a transparent 
separation of the different scales in the problem. As an 
application we have derived an exact all-order expression 
for the resummed deep-inelastic structure function F2, 
which is much simpler than corresponding results found 
in the literature. 
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